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Abstract 

We study spectral properties of a non-Hermitian Hamiltonian describing a 
quantum particle propagating in a random imaginary scalar potential. Cast in 
the form of an effective field theory, we obtain an analytical expression for the 
ensemble averaged one-particle Green function from which we obtain the den- 
sity of complex eigenvalues. Based on the connection between non-Hermitian 
quantum mechanics and the statistical mechanics of polymer chains, we de- 
termine the distribution function of a self-interacting polymer in dimensions 
d > 4. 
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The field of non-Hermitian quantum mechanics has, in recent years, attracted great 
interest. A variety of applications have been identified including the study of anomalous 
diffusion in random media |IJ, scattering in open quantum systems [0], neural networks ||, 
chiral symmetry breaking in quantum chromodynamics and the statistical mechanics 
of flux lines in superconductors ||. The last of these has involved the study of the quantum 
mechanics of a particle confined to a random impurity potential and subject to an imaginary 
vector potential. This work revealed a novel mechanism of "delocalisation" of the quantum 
particle, sharply contrasting with the behaviour of the Hermitian counterpart [|6|-|l^|. Here 



we investigate the spectral properties of a quantum particle confined to an imaginary scalar 
potential. To motivate our investigation, we apply these results to the study of the statistical 
mechanics of a self-interacting polymer chain. 

The Hamiltonian describing a particle propagating in a random scalar potential is defined 

by 

*=|U.V(r), (1) 

where the potential V is drawn from a Gaussian distribution with zero mean, and correlator 
(V(r)V(r'))v = 7<5 d (r — r'). To help motivate our discussion, we will apply our analysis to 
the statistical mechanics of a polymer chain. 

In the continuum limit, the probability distribution Z(r,t) of the end-to-end distance 
r of a self-interacting polymer chain of length t can be expressed in the form of a path 
integral [ff3Lff^, 



z ' r ' *> - Co Dx(r) exp {- II dT 'j (i) ' - £ II <^"I^m - • < 2 > 

The first term, the Wiener measure, determines the entropic contribution, while the second 
term represents the repulsive local contact interaction between the monomers that make 
up the chain. Decoupling the interaction by means of a Hubbard-Stratonovich field, the 
distribution function can be identified with the ensemble average of the Feynman propagator 
of the Hamiltonian above ([l|). Defining U(t) = exp[— tH], 
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the distribution is obtained from the ensemble average Z(r,t) = (U(r,t))y. 
Applying a spectral decomposition of the complex Green function 
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where \Ri) and (Lj| denote the right and left-hand eigenfunctions of H, and Zi denote the 
complex eigenvalues, the distribution function (0) takes the form 
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where the integration runs over the entire complex plane. Previous studies of the spec- 
tral properties of weakly non-Hermitian operators have largely (although not exclusively) 
focussed on properties of random matrix ensembles [p!^ -|T8|. Such studies [0] have em- 
phasised the pitfalls of a diagrammatic analysis based on a perturbative expansion of g in 
the random potential. The problems are revealed by representing the density of complex 
eigenvalues through the operator identity 
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wherein the Green function is shown to be non- analytic everywhere in which the density of 
states (DoS) is non-vanishing. To circumvent these difficulties, a representation has been 
introduced ]l7|,7fl in which the complex Green function is expressed through an auxiliary 
Hermitian operator, 

e z-H" 
z* - #t e 

where z = x + iy and cr represent Pauli matrices. Making use of this construction, a 
relationship between g and the matrix Green function is straightforwardly obtained, 
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g(z) = lim G 21 (e). 
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At the same time, this representation manifests an implicit time-reversal and chiral sym- 
metry of the matrix Hamiltonian: a 3 G^ 1 (0)a 3 = — G _1 (0). Recently, this representation 
was successfully combined with diagrammatic perturbation theory to study the spectrum 
of the Fokker-Planck operator describing particles diffusing in a quenched random velocity 
field [0. 

Cast in this form, we are able to apply standard field theoretic methods to obtain statis- 
tical properties of the single-particle Green function. Our approach is closely related to that 
recently developed by Efetov p(| (see also Ref. to investigate spectral properties of 



the random Schrodinger operator subject to an imaginary vector potential. This approach 
involves a generalisation of the supersymmetry method originally tailored to the description 
of disordered conductors ||22|| . Being somewhat more technical than the diagrammatic per- 



turbation theory employed in |[19|| , the use of the supersymmetry technique is nevertheless 
justified for the problem at hand. Indeed, although the diagrammatic approach reproduces 
the mean- field result, it does not take properly into account the existence of the massless 
Goldstone modes (see below), and therefore can not be safely used if one wants to go beyond 
the saddle-point approximation. 

The analysis begins by expressing the matrix Green function as a functional integral over 
8-component supervector fields G a p(0, r) = «(tr[i2 _1 \l/(0) ® 1 l r (r)cr3i?E /3a ])^/4, where 
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R = exp[— z7T<Ti/4], and the infinitesimal imaginary part ensures convergence. Here we 



adopt a standard notation [22 in which the fields \I/ subdivide into a time-reversal (tr), 
a Fermion-Boson (fb), and a "spinor" or matrix sector. S' 3 " is a 2 x 2 matrix which projects 
on to the a/3 components in the spinor space. 

Expressed in this form, an ensemble average over the random impurity potential gen- 
erates a quartic interaction of the fields which can be decoupled by the introduction of 
8x8 component supermatrix fields Q. Taking our notation from disordered conductors and 
defining t~ 1 (x) = 2^ r yv{x) ) where u(x) is the unperturbed DoS, we obtain 
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where strM = M FF — M BB represents the trace operation for supermatrices. The superma- 
trix fields Q have an algebraic structure which reflects that of the dyadic product \I/ Cg> 
Integrating over the superfields we obtain (G Qj g(0, r))y = — (tr[i? _1 ^(0, r)cr 3 _RS /3a ])Q/4 
with 
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where the supermatrix Green function is defined by 
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Further progress is possible only within a saddle-point approximation, which is controlled 
in the limit of weak disorder [y, 1/r] <C x. Minimising the action (|TT1) with respect to 
variations in Q, we obtain the saddle-point equation 



Q(r) = g(r,r) 



(13) 



The saddle-point solution is found from the ansatz that Q is homogeneous in space, and 
diagonal in the TR and FB sector. Applying the parametrisation Q = q + q ■ cr, we obtain 
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Qi = 12 = 0, g 3 = < 
^sgn(y) [ 0. 

Moreover, invariance of the Eq. (|1^) under rotations Q — > TQT^ 1 where [T, = 0, shows 

that the first solution of the saddle-point equation spans a degenerate manifold, Class CI 

in the classification of Ref. p4| . Below, we will find that this solution corresponds to the 

non-analytic part of the Green function with a non-zero DoS. The second (non-degenerate) 

solution yields the analytic part of the Green function, and does not contribute to the DoS. 

The region in the complex plane where the first solution is stable defines the support of the 

spectrum. Straightforward stability analysis shows that the boundary of the spectrum is 

determined by the equation y = l/r(x). 



Expanding the action in slow fluctuations T(r) around the saddle-point solution we 
obtain the low energy effective action 

S[Q] = ~ J drD(y) str (dQ)\ (15) 

where Q = Ta 3 T _1 , and D(y) = (1 — y 2 r 2 )D , with D = 2xr(x)/md, denotes the y- 
dependent classical diffusion constant. Thus, in contrast to a real random impurity potential, 
modes of density relaxation of the matrix Hamiltonian are controlled by massless Goldstone 
modes of a supersymmetric non-linear cr-model of symmetry class CI with a diffusion con- 
stant which depends explicitly on y. 

Expanding the supermatrix Green function, 

0(0, r) = -mv{x)f d {v) [q + q 3 Q(r)} , (16) 

where /d(r) = IrnGg (0, r)/ImGo (0, 0) and Gq = (x — p 2 /2m — i/2r)~ 1 , we obtain 

(G af3 (0, r)> = i™f d {r) J DQ str [iT 1 (q + Q(r)) a 3 R^ a ] e~ s ^. (17) 

Applied to the DoS, the projection matrix takes the form S = (o~i + io"2)/2. In this case the 
supermatrix degrees of freedom of Q are decoupled from the source, and the DoS is specified 
simply by the mean-field result 

f (47T7)" 1 , \y\ < 1/t(x), 
p{z) = (18) 
[0 \y\>l/r(x), 

satisfying the sum rule / dyp(z) = v(x). This result compares with that obtained for the 



corresponding random matrix ensemble |18|| , and contrasts that obtained for an imaginary 
vector potential |p0| . In particular, since the DoS source does not couple to the effective 
action, the mean field estimate is unchanged by integration over Q. As a result, the complex 
eigenvalue density remains non- singular at y — 0. 

We have computed numerically the eigenvalues of the two-dimensional lattice version of 
the Hamiltonian (j^). Theory and simulation are compared in fig. |T], for 50 realisations of 
a 32 x 32 lattice with 7 = 1. The solid line represents the boundary of the spectrum as 
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calculated in the saddle-point approximation, y = ±l/r(x) = 27rz/(x)7, —4 < x < 4, where 
v{x) is the exact DoS of the clean two-dimensional tight-binding model. The results show 
good agreement between theory and numerics, with deviations becoming significant only 
near the edges of the band, and in the vicinity of the band center. (Note that, in contrast 
to the Hermitian disordered Hamiltonian, the imaginary scalar potential raises the energy 
of the low lying states, an effect easily understood within the framework of second order 
perturbation theory. However, the localisation properties of such "Lifshitz tail" states, as 
well as their sensitivity to optimal fluctuations of the impurity potential is a subtle question 
which lies beyond the scope of the present investigation. Secondly, a sublattice symmetry 
specific to the tight-binding Hamiltonian induces a reflection symmetry of the spectrum 
around x = 0. This additional symmetry, which parallels that found with the imaginary 
vector potential, leads to the accumulation of a finite fraction of states along x = 0, a 
phenomenon which has no counterpart in the continuum model.) We believe that the small 
deviations of numerics and theory away from the band edge and band centre can be ascribed 
to the influence of massive fluctuations of the matrix fields, i.e. those which violate the 
symmetry [T(r),a 2 ] = 0. 

Turning to the self-interacting polymer chain discussed above, Eq. QIE| ) can be used 
to obtain the probability distribution Z(r,t). However, to do so, we should recall that 
the saddle-point approximation is justified in the limit xr ^> 1 which, setting m — 1, 
translates to the following condition on the contact interaction: '-yx^ -4 )/ 2 <C 1. The latter 
is satisfied as x — > in dimensions higher than four, and x — > oo in dimensions lower than 
four. Equivalently, applied to the polymer chain, the analysis above applies at time scales 
t > 7 2 /( d ~ 4 ) for d > 4, and t < (l/-f) 2 /^- d ^ for d < 4 defining the upper critical dimension 
as d c = 4 |TJj] . Evaluation of Eq. (|5]) for d > 4 gives the power spectrum 



Z( P ,t) = s MpM eM~P% w(p) = H P r. (19) 
uj{p)t 



This result parallels that obtained by Chalker and Wang [BJ for the time evolution 
of the particle density of the random Fokker-Planck operator. In particular, in the limit 



t — > oo it recovers diffusive behaviour (r 2 ) ~ t, while at short times, when the saddle-point 
approximation is uncontrolled, we have (r 2 ) ~ (^t) 2 ' ^ d ~ 2 ^ . The long-time behaviour for d < 4 
is strongly affected by low-energy states which, being sensitive to optimum fluctuations of 
the random impurity potential, is not accessible within the present framework. 

In conclusion, applying conventional field theoretic methods, we have shown that spec- 
tral properties of the non-Hermitian Hamiltonian describing a particle propagating in an 
imaginary scalar potential are governed by a supersymmetric non-linear cr-model (of sym- 
metry class CI ||24||). The latter has been used to obtain the complex eigenvalue density. 
As an application of these results, we have obtained the power spectrum of the probability 
distribution function of a self-interacting polymer chain. Although the analysis developed 
here has focussed on one-particle properties, an extension of the present approach to treat 
fluctuation phenomena is, at least formally, straightforward (see, for example, Ref. p5[). 
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FIGURES 

FIG. 1. Complex eigenvalues taken from 50 realisations of a 32 x 32 square lattice tight-binding 
model with 7 = 1.0. The solid line is y = ±l/r(x) = 2irv(x)'y. 
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